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Jh ■ Abstract 

It is shown that an effective anisotropic spherically symmetric fluid 
model with heat flow can absorb the addition to a perfect fluid of 
^- ■ pressure anisotropy, heat flow, bulk and shear viscosity, electric field 

and null fluid. In most cases the induction of effective heat flow can 
q ! be avoided. 



O 1 Introduction 

O' 

Spherically symmetric perfect fluid solutions in general relativity have been 
studied from its very beginning, starting with the interior Schwarzschild so- 
lution. Gradually different mechanisms in stellar models have been identified 
that create pressure anisotropy [I], [2J. Radiating spherical collapse demands 
the introduction of heat flow [3] or null fluid which describe energy dissipa- 
tion in different approximations [I]. More realistic fluids possess also bulk 
and shear viscosity [5], [6], [7]. Charged perfect fluids or dust have been 
discussed by many authors [8] . It has been shown too that the sum of two 
perfect fluids, two null fluids or a perfect and a null fluid can be represented 
by effective anisotropic fluid models [9]. In view of the many existing rela- 
tions among the fluid models mentioned above it is hard to point out some 
centre. 

In this paper it is shown that the anisotropic fluid with heat flow is in some 
sense the most fundamental model and can absorb the addition of viscosity, 



charge and null fluids. Something more, heat flow is not generated in most 
cases. 

In Sec. 2 the basic anisotropic fluid model is defined. In Sec. 3 the fluid 
is supplied with bulk and shear viscosity which leads to a new effective 
anisotropic model. In Sec. 4 the same is done for the addition of charge 
and in Sec. 5 null fluids are accommodated into the anisotropic model. Sec. 6 
summarizes all additions and the effective characteristics of the fundamental 
anisotropic model are given. Several conclusions are drawn. In Sec 7 the 
static case is discussed and a formula for the general solution is presented. 
Sec 8 contains some basic conclusions. 

2 Anisotropic fluid model with heat flow 

Einstein's field equations are given by 

8vrT a/9 = G a/3 (1) 

where G a p is the Einstein tensor, T Q/ g is the energy-momentum tensor (EMT) 
and units are used so that c — G — 1. The general spherically symmetric 
metric is written as 

ds 2 = -A 2 dt 2 + B 2 dr 2 + R 2 (d6 2 + sin 2 9 dip 2 ) (2) 

where A,B,R are positive functions of t and r only. The spherical coordinates 
are numbered as x° = t, x 1 = r, x 2 = 8, x 3 = (p. The Einstein tensor involves 
the Ricci tensor and scalar which are given by the metric and its first and 
second derivatives [4]. Its non-trivial components are Goo, Gm> Gn, G22 = 
sm- 2 9G 33 . 

We are interested in the structure of EMT. It reads for anisotropic fluids 
with heat flow 

T a p = (fj, + p t ) u a u,3 + p t g a(3 + (p r - p t ) XoXa + QaU/3 + u a q p . (3) 

Here [i is the energy density, p r is the radial pressure, pt is the tangential 
pressure, u a is the four velocity of the fluid (a timelike vector), \ a is a un h 
spacelike vector along the radial direction and q a is the heat flux (in the 
radial direction too). We have 

u a u a = -l, x a X a = l, u a Xa = 0, u a q a = 0. (4) 



It is assumed that the coordinates are comoving, hence, the fluid is motionless 
in them 

U a = A~ 1 ^, X a --B- X b^ q a = qB- x b1 (5) 

(6) 



where q = q(r,t). This gives 

Too = M 2 , Toi = -qAB, T n = p r B 2 , T 22 = p t R 2 , 
which should be plugged in the Einstein equations. They become 
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Here the dot means time derivative and the prime is a radial derivative. 

The anisotropic fluid does not radiate when q = and becomes perfect 
when p r — pt- Thus it accommodates anisotropy of pressure and heat flow 
when one starts with a perfect fluid model. 

Let us see now what happens when other EMTs are added to the basic 
anisotropic one. 



3 Bulk and shear viscosity 

Bulk viscosity j5], [6], [7] adds to the basic EMT the following piece 



(11) 



where ( is a coefficient, G is the expansion of the fluid and h a p is the projector 
on the hyperplane orthogonal to u a 

6 = u°t a , hap = g a p + U a Up. (12) 

Obviously, this means the appearance of effective pressures 

Vr = Pt = -C©, (13) 

which should be added to p r ,Pt- They do not change the degree of anisotropy 
Ap = p T — Pt- Thus even perfect fluid can absorb the bulk viscosity The 
quantities fi, q remain the same. No heat flow is generated in particular. 
Shear viscosity is responsible for the piece 

T% = -2 V h ai h^ 8 (14) 

where <7 Q/ g is the shearing tensor 

0"ay8 = M (a;/3) + a (a M /3) - 77®^' (^) 

i] is some coefficient and a a is the acceleration 

a a = u a tfw . (16) 

The shearing tensor satisfies the conditions 

a^u 13 = 0, a a/3 g^ = 0, (17) 

hence, Eq (14) transforms into 

T^ = -2 WaP . (18) 

Use of Eqs (5,15) gives the non-zero components of the shear [4] 

m = |aV, a 22 = ^¥- = -^V (19) 

6 sin 6 

where 



\o 2 = a^aap. (20) 

One can check that the same components follow when a a is written as the 

tensor 

1 

<Ta$ = -^h a/3 + aXaX/3- ( 21 ) 



Thus it coincides with the general shear tensor defined by Eq (15) in the 
spherically symmetric case. It also satisfies relations (17) in any metric. 

Plugging Eq (21) into Eq (18) and comparing it to Eq (3) we find the 
effectively generated pressures 

p s r = -2pf = ~na. (22) 

The degree of anisotropy is changed. There is no generation of energy density 
or heat flow. The scalars G, o can be expressed through the metric and its 
first derivatives: 

^ 1 ( B k\ 1 ( B h\ 

6 = a[b +2 r)' "Ab-rY (23) 



4 Electromagnetic fields 

The EMT of electromagnetic fields is given by 

T^ M = ^ (f^ - \g a ,F, u F^ (24) 

where F^ v is the Faraday tensor. One defines a unit timelike vector field n M . 
An observer moving in its direction will measure electric and magnetic field 
respectively 

These fields are spacelike, E a n a = H a n a = 0. The Faraday tensor decom- 
poses like [TO] 

F a /3 = ZaPnH^ — 2_E'[ Q n / g] . (26) 

Plugging this expression into Eq (24) gives formula (7) from Ref |10J . which 
becomes after some rearrangements 

T™ = ^{H 2 + E 2 ) (n a n p + ± 9o ^\ - i- (E a E p + H a H fi )+2j (a n fi) (27) 

where E 2 = E^E 11 , H 2 = H^H^ and j a is the Poynting vector that measures 
the energy flow in the spacetime 

Ja = -U^^/T. (28) 

47T 



In order to absorb this EMT by the EMT for anisotropic fluid we choose 
the direction n a = u a and x a = E a /E. The latter is possible because when 
spherical symmetry is imposed H a = and E a has only a radial spatial 
component. The would be heat flow term in Eq (27) disappears and we get 

E 2 
T ap = 2eu a u p + eg aP - 2ex a X/3, e = ^~- (29) 

Thus the addition of electric field induces effective pressures and energy den- 
sity 

/i £ = pf = ~Pr = e, (30) 

related by simple linear equations of state. Hence, a charged perfect fluid or 
a charged anisotropic fluid may be represented effectively by some neutral 
anisotropic fluid. There is no heat flow induction in this case. 

Charged fluids, however, must satisfy the Maxwell equations in addition 
to the Einstein ones. It is clear from the above that the electromagnetic 
tensor has a single component. The first pair of Maxwell equations gives 
F w = $', where $ is the only component of the electromagnetic potential. 
The second pair of equations reads 



,a/3 = /_„\-l/2 



1/2 ™/3 



4nTu a = F ap . p = (-g)~ 1/z (- g y /z F ap (31) 

where r is the charge density, g is the determinant of the metric, usual deriva- 
tive is denoted by comma and the covariant derivative is denoted by semi- 
colon. In the spherically symmetric case this formula provides two equations. 
One of them gives 

|^ = P (r), (32) 

P (r) being an arbitrary function of the radius. Hence, the combination in 
the l.h.s. depends only on the radial coordinate. Then the second equation 
yields 

4- = IL (33) 

and becomes a formula for the charge density. Plugged in the definition of 
E the Maxwell equations lead to 

P(r)' 
E = ^ (34) 
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Substituting Eq (36) into 


Eq (35) one 


finds 





which represents a constraint on the form of the effective density and pres- 
sures e in the general time-dependent case and no constraint in the static 
case. 

5 Null fluid 

Null fluid describes dissipation in the free streaming approximation jl] and 
adds to the basic EMT the piece 

(35) 
(36) 
(37) 



T a/3 = £U * U P + £ XaX(3 + £ {UaXp + UpXa) ■ ( 38 ) 

A comparison of this expression with Eq (3), taking into account that q a = 
qx a , shows that the addition of null fluid generates effective energy density, 
radial pressure and heat flow, all of them equal 

A*" = Pr =q N = e. (39) 

No tangential pressure is generated. This is the only case where an effective 
heat flow is induced. 



6 Summary 

The results in the previous sections show that viscosity, electric charge and 
null fluids are equivalent to induced energy density and pressures, related 
by simple linear equations of state \x = np r , n = 0, ±1 and p t = kp r , k = 
0, ±1,-1/2. When all such additions are combined and absorbed by the 
initial anisotropic fluid model, one obtains an effective model with 

// = fM + e + e, (40) 



4 
Pr = Vr - C© - 3 W - e + e, (41) 

p* = ^ - (6 + -rja + e + e, (42) 

q e = q + e. (43) 

These should be plugged into Eq (6) and hence in the l.h.s. of the Einstein 
equations (7-10). There are 4 equations for 8 functions (fj,,p r ,p t ,qX,Vi e i £ )- 
The quantities a, are given by Eq (23). One has to impose 4 relations on 
these functions or set some of them to zero in order to obtain a determined 
system of equations. Several conclusions can be drawn. 

Only the functions giving the two modes of dissipation of energy (q, e) 
have effect upon the heat flow. 

Viscosity ((, rj) does not induce effective energy density. 

An important characteristic of the fluid model is the anisotropy factor 

A p e = p e r - p e t = Ap - 2r](T - 2e. (44) 

We see that shear viscosity and charge induce pressure anisotropy. Their 
absorption by a perfect fluid (Ap = 0) makes the latter anisotropic and adds 
several more sources of anisotropy to the usual ones [lj. 

Bulk viscosity and null fluid don't induce anisotropy and may be absorbed 
by an effective perfect fluid model with heat flow. 

Charged perfect fluids are equivalent to anisotropic neutral fluids because 

Ap=-2e. (45) 

In addition to the Einstein equations charged fluids satisfy Maxwell equations 
but we have clarified that they give a formula for the charge density and a 
mild constraint on e. Thus results about anisotropic fluids may be carried 
over to charged perfect fluids and vice versa. 



7 The static case 

In this case there is no time dependence and R = r, q = 0. The other three 
Einstein equations reduce to 
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The difference of the last two equations gives a linear equation for B 2 . 
Introducing the variable z 

A 2 _ e 2f(z-2/r)dr ^ 

and solving the linear equation one obtains [UJ 

2 z 2 K 2 , , 

B 2 — (50) 

r 6 [2 / jg (1 + 8vrApr 2 ) K 2 dr + C] 
where C is a constant of integration and 

K = e ^^ +z ) dT . (51) 

Thus the two generating functions Ap and z determine the metric of any 
anisotropic fluid solution. Eqs (46-48) then determine its energy density and 
the two pressures. In particular, using the results in the previous sections, 
we can plug formulas (40-42) into the left hand sides of Eqs (46-48) to obtain 
the general solution of any fluid with anisotropy, shear and bulk viscosity 
and electric charge. 

8 Conclusions 

Anisotropic fluid models were studied for the first time by Lemaitre in 1933 
on Einstein's recommendation (see the Golden oldie [12] and the references 
therein). His work, however, had no impact for a long time and perfect 
fluid models have been examined for decades. In later times physical reasons 
were given for anisotropy of pressures in relativistic star models PQ, [3] and 
now there exists an extensive literature on this topic. In the present paper 
we have given additional arguments in favour of the anisotropic fluid model 
with heat flow as an effective model, encompassing fluids with anisotropy, 
viscosity, charge and radiation. We have also given a vocabulary of the 
characteristics of different fluid models in terms of the basic anisotropic one. 



Hence, one can concentrate on finding its solutions and then translating them 
to the other models. This procedure should yield numerous, albeit formal, 
relations between the solutions of anisotropic and charged perfect models, or 
between viscous and anisotropic models and so on. Such direction of research 
deserves a lot of further study. 
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